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We calculate the properties of the Ad ferromagnet SrRuC>3 in bulk and thin film form with the aim 
of understanding the experimentally observed metal to insulator transition at reduced thickness. Al- 
though the spatial extent of the Ad orbitals is quite large, many experimental results have suggested 
that electron-electron correlations play an important role in determining this material's electronic 
structure. In order to investigate the importance of correlation, we use two approaches which go be- 
yond the conventional local density approximation to density functional theory (DFT): the local spin 
density approximation + Hubbard U (LSDA+C) and the pseudopotential self-interaction correction 
(pseudo-SIC) methods. We find that the details of the electronic structure predicted with the LSDA 
do not agree with the experimental spectroscopic data for bulk and thin film SrRuOa. Improve- 
ment is found by including electron-electron correlations, and we suggest that bulk orthorhombic 
SrRuC>3 is a weakly strongly- correlated ferromagnet whose electronic structure is best described by a 
0.6 eV on-site Hubbard term, or equivalently with corrections for the self-interaction error. We also 
perform ab initio transport calculations that confirm that SrRuOs has a negative spin polarization 
at the Fermi level, due to the position of the minority Ru Ad band center. Even with correlations 
included in our calculations we are unable to reproduce the experimentally observed metal-insulator 
transition, suggesting that the electronic behavior of SrRuOa ultra-thin films might be dominated 
by extrinsic factors such as surface disorder and defects. 

PACS numbers: 71.20.-b, 79.60.Dp, 79.60.-i, 32.10.Dk 



I. INTRODUCTION 

Conductive oxides are essential components in compos- 
ite oxide heterostructures where they are of ten u sed as 
electrode materials in thin film applications .1112121 l n the 
perovskite crystal family (ABO3 stoichiometry) , the itin- 
erant ferromagnetic S1RUO3 (SRO) is a popular choice 
since it is one of the more conductive metallic oxides 
with good thermal properties!^ In thin films, SRO is in- 
tensely investigated as a possible route to the realization 
of novel field-effect devices In addition, it is of par- 
ticular interest to the spintronicP^ and multiferroicP^ 
communities, which have been recently energized by the 
possible device applications available from engineering 
interface phenomena! 11 ! 12 ! 13 ! 14 ! 15 ! 16 -^ However, one lim- 
itation in the design of thin film oxide devices is the ob- 
servation of increased resistivity in metal oxides as the 
film thickness decreases. Such behavior is clearly present 
in ultra-thin films of SrRu03, where a metal-to-insulator 
(MI) transition!^ occurs at four monolayers. This sub- 
stantial change in the electrical conductivity presents a 
serious challenge for device miniaturization. In this work 
we explore the underlying physics of the thin film MI- 
transition, which to date remain to be understood. 

The 3d transition metal oxides (TMOs) are known 
to possess strong electron-electron correlation effects 
that can drive a system that should be metallic 
within a simple band picture into an insulating state. 
Due to the large spatial extent of the 4<i-orbitals 
in the ruthenates, correlation effects are anticipated 
to be less important as stronger hybridization pro- 



vides more effective screening and a reduced Hubbard 
U (Coulomb repulsion energy). Many experimental 
studies have already addressed the degree of electron- 
electron correlation in SrRuOa including X-ray and ul- 
traviolet photoemission spectroscopy! 1 ^- 9 ! 20 ! 21 ! specific 
heat measurements^^ infrared and optical conductiv- 
ity measurements,^ and transport experiments.^ For 
example, Kim and coworkers^! use X-ray photoemis- 
sion spectroscopy (XPS) to identify how such correla- 
tions change within the ruthenate family, and Toyota 
et al^ use photemission spectroscopy (PES) to detail 
the metal-insulator transition in SrRuOa as a function 
of film thickness concomitant with the onset of mag- 
netism. In all of these studies, the general consensus 
is that electron correlation effects do play a role in de- 
termining the electronic structure of this itinerant ferro- 
magnet, but to what degree remains unclear. Further- 
more, some theoretical investigations have begun exam- 
ining covalencyj^S correlation^ and orbital ordering^ ef- 
fects in bulk SrRuOa. The magnetic properties of SRO 
under epitaxial strain have also been investigated with 
first-principles techniques! 2 ^! 

In this work, first-principles density functional theory 
(DFT) calculations are performed, first to identify the 
degree of correlation in bulk SRO, and second to inves- 
tigate the driving force for the metal-insulator transition 
in ultra-thin films. We use two approaches to introduce 
correlation effects into the conventional band theory (lo- 
cal spin density) approach for treating the Ru 4d-orbitals 
and their hybridization with O 2p-orbitals: the local 
spin density + Hubbard U (LSDA+f/), and the pseu- 
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dopotential self-interaction corrected (pseudo-SIC) local 
spin density methods. In addition we investigate two 
structural variants - the ideal cubic perovskite structure 
and the experimentally observed orthorhombic structure, 
which includes tiltings and rotations of the RuC^ octa- 
hedra. By comparing the results that we obtained for 
both methods and both structure types, we are able to 
comment on the nature of the metal-insulator transition 
in ultra-thin films. 



II. 



CRYSTAL STRUCTURE & MAGNETISM 



The perovskite class of materials is described by a net- 
work of corner-sharing BOq octahedra, in which the A- 
site cation is located at the center of a cube defined by 
eight BOq units. The ideal perovskite is cubic (space 
group Pm3m) , however several modifications exist owing 
to the range of cation sizes that can be accommodated in 
the structure. Deviations from the ideal cubic structure 
are defined by the Goldschmidt tolerance factor 



t' 



Ra + Ro 
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where Ri is the radius of atom i, and can be attributed to 
the requirement to optimize the anion coordination about 
the A-site cationPS Using the Shannon-Prewitt radii for 
this compound, a predicted tolerance factor of £'=0.908 
is found, which is far from the ideal case t'=l, suggesting 
that distortions should occur. 

Indeed, SRO undergoes a series of structural trans- 
formations with temperature, from high symmetry cubic 
(Pm3m, stable above 950 K) to tetragonal (I4/mcm, sta- 
ble between 820 K and 950 K) to distorted orthorhombic 
structure (Pbnm) at low temperatures. The orthorhom- 
bic distortion from the ideal cubic can be described by 
the tilting of the RuOe octahedra in alternate directions 
away from the c axis, and the rotation of the octahe- 
dra around the b axis; in both cases adjacent octahedra 
distort in the opposite sense (Figure fl|. The degree of 
tilting and rotation (as defined in Fig. [l| of the octahe- 
dra are useful in describing the distortions in the oxygen 
network from the perfect cubic case. A rotation angle 
of 7.56° and a tilting angle of 10.47° (corresponding to 
a Ru-O-Ru angle of 159°) are found for Pbnm SrRuC>3. 
The structural changes reduce the hybridization between 
the Ru Ad states and O 2p states and lead to a narrowing 



of the bandwidths (see Section IV I compared with the 



ideal cubic case. Consequently the degree of correlation, 
described by U/W where W is the valence bandwidth, is 
expected to be enhanced. 

Below approximately 160 K, SrRuC>3 exhibits strong 
ferromagnetic behavior, and has a measured Rhodcs- 
Wohlfarth raticpSl (fi e g/fi sat ) of 1.3 suggesting that its 
magnetism can be well described by a localized d-electron 
model similar to the elemental ferromagnetic metals. 
Within this model and under an octahedral crystal field, 




FIG. 1: (Color online) The orthorhombic (Pbnm) crystal 
structure of SrRuOa. The unit cell contains four formula units 
(f.u.) of the ideal cubic perovskite (Pm3m). The structure 
is stabilized by shortening the Sr-O distance followed by a 
cooperative distortation of the RuOe octahedra to reduce the 
coordination volume of the Sr ions; this results in a smaller 
Ru-O-Ru bond angle, which in turn decreases the Ru 4d 
bandwidth and metallicity. The deviation in the structure 
from the high symmetry cubic state can be quantified using 
the tilting angle (18O°-0)/2 and the rotation angle (9Q°-0)/2 
of the oxygen octahedra. 



the 4c?-manifold splits into a threefold degenerate t2 g sub- 
band that is lower in energy than the twofold degener- 
ate e g band. Neglecting covalency, we would expect a 
spin-only magnetic moment of 2 [Ib , corresponding to a 
low-spin state for the Ru 4+ ions (<i 4 : t\ g ^,t\ g ^}. Experi- 
mentally, however, the moment is measured to be closer 
to 1.1 /is/f.u., although values ranging from 0.9 /is/f.u. 
and 1.6 /Us/f.u. have also been reported.^ (The spread 
in values is attributed to the large magnetocrystalline 
anisotropy of the material, and the difficulty in mak- 
ing large single-domain samples.) First-principles cal- 
culations also report a ma gnetic mo ment ranging from 
0.9 [x B /f.u. to 2.0 Ms /f. u piMl31H The reduced calcu- 
lated magnetic moment in the solid compared to that in 
the free ion limit is due in part to the large spatial extent 
of the Ru 4c? orbitals, which results in a significant over- 
lap (hybridization) with the oxygen 2p. Furthermore, 
due to the metallic character of SRO, an overlap of the 
majority and minority Ru 4d bands occurs at the Fermi 
level; as a result partial occupation of the minority band 
also leads to a reduced magnetic moment. 

In this work, we examine the LSDA and "beyond- 
LSDA" electronic and magnetic properties of both the 
Pbnm and Pm2>m crystal variants. Metallicity and mag- 
netism are both related to the d-band width, which in 
turn depends on both correlations and structural prop- 
erties such as tiltings and rotations of the oxygen oc- 
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tahdera. Our goal, therefore, is to identify the relative 
contributions of electron-electron correlation effects and 
structural distortions in driving the metal-insulator tran- 
sition in SrRu03 thin films. 



III. THEORETICAL METHODS 
A. LSDA 

Our initial electronic band structure calculations were 
performed within the local spi n densit y appro ximati on^ 
(LSDA) using both the siestaI^SMI^H and VASr^Uini den- 
sity functional theory (DFT) packages. In each method 
we used the Perdew-Zunger-^Q parametization of the 
Ceperley-Alder data 42 for the exchange and correlation 
(XC) functional. 

The core and valence electrons were treated with the 
projector-augmented-wave (PAW) methooP for all cal- 
culations performed with VASPpS Furthermore a plane- 
wave energy cutoff of 500 eV was used and found to 
produce excellent convergence. The three-dimensional 
Brillioun zone was sampled with a 12 x 12 x 12 fc-point 
Monkhorst-Pack mesrfHl for the cubic bulk (5 atom unit 
cell) and thin film structures and al2xl2xl0 fc-point 
sampling mesh for the bulk orthorhombic structure (20 
atom unit cell). For the orthorhombic films we used a 
12 x 12 x 3 fc-point sampling. In all cases the tetrahe- 
dron method with Blochl correctional was used for the 
Brillouin zone integrations. 

In the localized basis code SIESTA, the core and va- 
lence electrons were treated with norm-conserving fully 
separable^ Troullier-MartirPS pseudopotentials.^ The 
localized atomic orbitals for each atom used a single-C 
basis set for the semicore states and a double-C for the 
valence states. Total energies were computed on a uni- 
form real space grid with a cutoff of 800 Ry in order 
to reach comparable accuracy to the planewave codeP^I 
The Brillioun zone of the cubic structure was sampled 
with a 26 x 26 x 26 fc-point Monkhorst-Pack mesh, while 
the Pbnm structure was sampled with a 15 x 15 x 12 
fc-point mesh. Integrations were performed with a Gaus- 
sian broadening of 0.10 eV in all calculations. 

The equilibrium lattice parameter for the cubic struc- 
ture was found by fitting the total energy as a function 
of volume to the Murnaghan equation of state. Excellent 
agreement was found between the two codes (Table [I]) , 
with a slight underestimate of the experimental lattice 
constant typical for the LSDA. The cell parameters and 
atomic positions of the orthorhombic structure (Table II]) 
were optimized by starting from the positions reported in 
Ref. [55] and the ionic coordinates were relaxed until the 
Hellmann-Feynman forces on the atoms were less than 
4 meV A -1 . 

Investigations of the electronic structure of cu- 
bic S rRu03 have been described by several different 
groups^* 28 * 32 * within the LSDA; here we briefly summa- 
rize their conclusions and remark that our results are con- 



a/a a 
B (GPa) 
B' 

Moment (/ib/f.u. 



VASP 

0.98 
200 
4.6 
1.09 



SIESTA 

0.98 
219 
4.4 
1.26 



TABLE I: Results obtained for cubic SrRuOs within the lo- 
cal spin density approximation (LSDA) from the two codes 
used in this work. The equilibrium lattice constant relative 
to the experimental value (ao = 3.9735 A) determined from 
high-temperature neutron diffraction dataj^ the bulk modu- 
lus B, the pressure derivative B' , and the magnetic moment 
per formula unit are given for each code. Very good agree- 
ment is found between the two codes. The notation for the 
codes are as follows, SIESTA: Spanish Initiative for Electronic 
Simulations with Thousands of Atoms local orbital code and 
VASP: Vienna Ab-initio Simulation Package planewave code. 



Atom 


Site 


X 


y 


z 


Sr 


4c 


-0.0050 


0.03039 


0.25 


Ru 


4a 


0.5 


0.0 


0.0 


0(1) 


4c 


0.0650 


0.4942 


0.25 


0(2) 


8d 


0.7158 


0.2834 


0.0340 



TABLE II: Calculated structural parameters for SrRuOs with 
the Pbnm symmetry using VASP. Our calculated orthorhom- 
bic lattice constants are a — 5.4924, b — 5.4887, and c = 
7.7561 A. 



sistent with the earlier calculations. A complete compari- 
son of the electronic structure of cubic SrRuOa calculated 
with both VASP and SIESTA is made in Section [IVJ In all 
cases, a metallic ferromagnetic ground state is found to 
be stable, with strong Ru 4d character at the Fermi level. 
Substantial hybridization occurs between the O 2p states 
and the Ru 4d-states and no energy gaps are observed in 
the densities of states. The calculated magnetic moment 
is also always reduced from the fully ionic limit of 2 /i^ . 



B. LSDA+(7 

The first "beyond-LSDA" method that we use to 
treat the exchange and correlation (XC) within DFT 
is the local spin density approximation with Hubbard 
U (LSDA+/7).H21 Here we use the spherically averaged 
form of the rotationally invariant LSDA+J7 introduced 
by Dudarev et aZ.J^in which only one effective Hubbard 
parameter, U c k = U — J, is used, where U and J are 
the spherically averaged Hubbard repulsion and intra- 
atomic exchange for electrons with the angular momen- 
tum of interest; in this case Ru 4<i states. We treat the 
double-counting term within the fully-localizcd-limit and 
note that this should more correctly describe the insulat- 
ing side of the metal-insulator transition that we study 
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here; an improved description of the metallic side might 
be achieved using the recently introduced interpolation 
between the around-mean-field and fully-localizcd-limit 
extremes.^ 

Within these approximations, the LSDA+C/ correction 
to the LSDA potential is 

A^(mm'a) = -(U - J) (p^, - \s m1 J\ , (1) 

where m, m! are the orbital indices, a is the spin index, 
and p a mm , is the orbital occupation matrix. The effect 
of the LSDA+C7 correction given by Eq. [T] is particularly 
transparent in the limit of diagonal p a mm i with orbital 
occupancies 1 or 0: Occupied orbitals experience a po- 
tential which is lower in energy by (U — J)/2 compared 
with the LSDA, and the potential for unoccupied orbitals 
is raised by the same amount. 

In this study we varied t/ e ff from to 6 eV for 
the Ru d-states, (the standard LSDA corresponds to a 
U c ff — eV). Structural optimizations were also per- 
formed within LSDA+C/ approximation, however negli- 
gible structural changes compared with the LSDA were 
observed. 



C. Self Interaction Corrections 

Our second approach for extending the treatment of 
the exchange and correlation is to correct for the spu- 
rious self- Coulomb and self-exchange interactions which 
arise within the LSDA, using the pseudo potentia l self- 
interaction corrected (pseudo-SIC) method!- 52 * 53 * 54 * These 
self-interaction errors are small for materials with delo- 
calizcd electronic states, but can be significant in systems 
with localized electrons where the interaction of an elec- 
tron with itself is large. 

Since the SIC in a periodic, extended system is not 
uniquely defined, many different methods have been pro- 
posed to remove SI in DFT calculations for solids (for 
a review see Ref. |55|). The procedure followed in the 
pseudo-SIC method is to: 

1. Project the occupied Bloch states onto the basis of 
the pseudo-atomic orbitals. 

2. Correct the potential for each Bloch state by the 
SIC for the pseudo-atomic orbital weighted by the 
projection, and scaled to account for the relaxation 
energy. 

Note that only the valence bands are corrected since the 
empty conduction bands, derived from orbitals where 
the occupation numbers are close to zero, are not self- 
interacting. This is in contrast to the LSDA+[7 method, 
in which the occupied bands are lowered in energy and 
the unoccupied bands raised. In principle, however, the 
two formalisms would yield equivalent results if a suit- 
able U — J (corresponding to the SIC energy) were ap- 
plied to all orbitals in the LSDA+C/ calculation. Indeed, 



whether the deficiencies of LSDA for strongly correlated 
systems derive from the absence of Hubbard U or the self- 
interaction error or both remains an open question. The 
pseudo-SIC method has some advantages over LSDA+J7, 
since it does not require a choice of which orbital to cor- 
rect, nor of the U or J parameters, and it can be applied 
readily to both magnetic and non-magnetic systems. We 
note that this is the first application of pseudo-SIC to 
an itinerant-correlated system, so the comparison with 
our LSDA+C/ results provides a test for the pseudo-SIC 
method. 



IV. RESULTS & DISCUSSIONS 

As we have mentioned, the electronic structure of 

SrRuOahas been investigated previously using the 
LSDA |22i28i23 Here we firgt revigit the LgDA with Qur 

own calculations, with an emphasis on understanding dis- 
crepancies bet ween the experimental measured photoc- 
mission result a 18 * 19 ! and the calculated LSDA electronic 
structure of the orthorhombic material. We then ex- 
tend our study to the two "beyond LSDA" approaches 
to examine correlation effects in bulk Pbnm and Pm3m 
SrRuOs. Finally, we examine unsupported films of 
SrRu03 in both structures in order to analyse the nature 
of the experimentally observed metal-insulator phase 
transformation . 

Cubic LSDA The total energies were calculated for 
both the ferromagnetically ordered and non-magnetic 
states of cubic SrRuOa using the optimized lattice pa- 
rameters. With both electronic structure codes the fer- 
romagnetic groundstate is always found to be lower in 
energy (vasp: 11.5 meV, SIESTA: 31.6 meV). Both codes 
yield very similar electronic structures. The density of 
states obtained using the vasp code is shown in Figure 
[2] The valence band is composed largely of O 2p states 
hybridized with Ru Ad states, with oxygen states predom- 
inately found in lower regions of the valence band and Ru 
states dominating at the Fermi energy. The large peak 
in the DOS near the Fermi level is caused by the fairly 
flat Ru t2 g bands near the Fermi level while the strongly 
hybridized e g orbitals form broader bands at the bottom 
of the valence and conduction bands. The Sr Ad states 
are found around 5 eV above the Fermi energy. 

The exchange splitting causes an energy shift between 
the majority spin and minority spin states; at the T- 
point a splitting of «0.50 eV is observed in the Ru Ad 
states, and of «0.20 eV in the O 2p states. The calcu- 
lated magnetic moments per formula unit are found to be 
1.09 pb (vasp) and 1.26 ps (siesta). With both imple- 
mentations of DFT, approximately 70% of the moment is 
found on the Ru atoms, with the remaining distributed 
about the oxygen network. The slight enhancement of 
the magnetic moment calculated with siesta compared 
to VASP is due to siesta's slight downward shift in en- 
ergy of the Ru t 2g band relative to the Fermi energy. 
We note that the consistency between the two DFT fla- 
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FIG. 2: (Color online) The total (gray line) and partial 
(shaded) spin-resolved densities of states for cubic SrRuC>3 
calculated within the LSDA using VASP. (upper) Sr id states, 



FIG. 3: (Color online) The total (grey line) and partial spin- 
resolved (shaded) densities of states for orthorhombic SrRuOa 
calculated within the LSDA using VASP. (upper) Sr 4d states, 



(middle) Ru 4d states [the t 2g and e g (unshaded, bold line) (middle) Ru 4d states [the t 2g and e 9 (unshaded, bold line) 



symmetries are shown], and (lower) O 2p states. The dashed 
line at eV denotes the Fermi level. 



symmetries are shown], and (lower) O 2p states. 



vors is essential to our later discussion of the effect of 
electron-electron correlations in the electronic structure 
of SrRuC>3, since the LSDA+J7 approach has been imple- 
mented in the VASP code, and the pseudo-SIC method in 
the siesta code. 

Orthorhombic LSDA Using the optimized LSDA lat- 
tice parameters for the Pbnm structure, we find that the 
ferromagnetic ground state is 6.34 meV/f.u. lower in en- 
ergy than the constrained paramagnetic structure, and 
additionally is 188 meV/f.u. (vasp) and 150 meV/f.u. 
(siesta) lower in energy than the ferromagnetic cubic 
phase. This energy stabilization can be associated with 
the oxygen octahedral tiltings and rotations, and agrees 
well with previous first-principles studieiP^^l that used 
experimental lattice parameters^ (the LSDA underesti- 
mates the lattice parameters by only about 1%). Also it 
is consistent with the experimental observation of ferro- 
magnetic SrRuC>3 in the distorted GdFeC"3 structure.^ 

The (P)DOS for the orthorhombic structure are shown 
in Fig. |3j and can be seen to be similar to those of the cu- 
bic structure discussed earlier (Fig. Consistent with 
the reduction in Ru Ad - O 2p overlap resulting from the 
tiltings and rotations, the bandwidths are slightly nar- 
rower in the orthorhombic structure, with the t^g band- 
width reduced by 0.35 eV, the e g by 1.5 eV and the O 2p 
by 0.60 eV. This results in a pseudo-gap opening in the 
minority e g bands at «-2.2 eV and a 0.20 eV gap open- 
ing «0.80 eV above the Fermi level. Interestingly, the Ru 
id exchange splitting is reduced slightly to 0.30 eV at T. 
This is accompanied by a reduction in the magnetic mo- 
ment compared with the cubic structure, to 0.79 /zs/f.u. 
(vasp) or 0.92 /ig/f.u. (siesta). As noted, the sponta- 



neous magnetization in the bulk (films) has been reported 
to be near 1.6 fis (1-4/xb), the LSDA underestimate is 
likely the result of the usual LSDA overbinding leading 
to enhanced Ru id - O 2p covalency. Note that this un- 
derestimate as compared to experiment!^ is significantly 
larger than that usually found for 3d transition metal 
oxides. It is further notable, because the orbital angular 
momentum is expected to be strongly quenched for Ad 
orbitals due to the cubic crystal field. We comment on 
the effect of including spin-orbit coupling later. 

Finally for the LSDA section, we compare our first- 
principles LSDA res ults with recent photoemission spec- 
troscopy (PES) dataP^EIMl w ith the goal of identifying 
which features are driven by correlation. In an ideal 
single-electron picture, the measured PES would consist 
of narrow peaks corresponding to the energies required 
to excite non-interacting electrons from the valence band 
into the continuum. However, the photoemission energies 
are more accurately interpreted as differences between 
two many-body A^-electron states: the ground state, and 
the excited state with a photoclcctron and hole. The 
effect of the many-body interactions is to broaden the 
one-electron peaks and shift spectral weight into so-called 
quasiparticle peaks. The strongest reduction in spectral 
weight from correlation effects occurs from so-called co- 
herent peaks near €f , and is accompanied by transfer 
of the spectral weight to higher energy features (inco- 
herent peaks) . Redistribution of the incoherent spectral 
weight into a well-defined satellite structure is indicative 
of strong correlations, whereas a redistribution into the 
background spectral distribution with a renormalization 
of the bandwidth indicates weak correlations. 



6 



i ' i ' i ' i 




• poly-crystal (hv = 40.8 eV) 
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FIG. 4: (Color online) The experimental PES spectra for 
bulk polycrystallineP 1 SrRu0 3 (filled circles) and 100 mono- 
layer SRO filrrP (triangles) grown on SrTiC>3 are compared to 
the calculated LSDA(+C/) and pseudo-SIC densities of states. 
The calculated DOS are broadened with an energy depen- 
dent Lorentzian (FWHM = O.f \e — e^j eV) and a Gaussian 
function (0.34 eV FWHM). An energy dependent parabolic 
background has also been added. 



In Figure [4] we show two experimentally measured 
spectra (see Ref s . H"8l and l59l for further information on the 
sample preparation) and our calculated LSDA results. 
First we comment on the discrepancies between the bulk 
polycrystalline spectrum and that of the thin film. Com- 
paring the experimental spectra, we sec that the thin film 
at 100 monolayers (which is representative of the bulk 
material) shows stronger coherent peaks than the highly 
broadened structure of the polycrystalline sample. In 
both cases however, the spectra are dominated by three 
principle features from the Ru ti g states near the Fermi 
level and the O 2p states between -8 and -2 eV. If we look 
more closely around the Fermi level, the polycrystal sam- 
ple has substantially reduced spectral weight, whereas 
near -1.3 eV it is enhanced compared to the film. This 
shift in the spectral weight to the incoherent peak agrees 
well with previous experimental comparisons^ made be- 
tween SrRuOs films and polycrystals and is attributed to 
the creation of near surface states induced during in situ 
scraping and not due to intrinsic correlation effects. Ad- 
ditionally the presence of grain boundaries and compo- 
sitional defects are also known to yield reduced coherent 
peak features in polycrystalline samples. For the remain- 
der of this study, we therefore restrict our comparison of 



the PES data to the 100 monolayer film, since it more 
accurately describes the intrinsic electronic structure. 

By comparing to our calculated densities of states, we 
can assign the peak features in the PES data to the cor- 
responding electronic states. In order to make the com- 
parison with our calculated DOS, we convolute an energy 
dependent Lorentzian function [full width at half maxi- 
mum (FWHM) = 0.1 je — ejr| eV)] with the calculated 
DOS to account for lifetime broadening. A Gaussian 
function with a FWHM — 0.34 eV is also used to account 
for the instrumental resolution, and an energy dependent 
parabolic background is also added. We find that the 
bands between approximately -2.5 eV and -8 eV are due 
to the O 2p states, with the peak of the non-bonding 
state centered at -3 eV. In the range between -8 eV up 
to the Fermi level are the occupied Ru 4c? states, with 
the tig state lying across the Fermi energy beginning at 
-3 eV. 

Some important discrepancies exist between our LSDA 
results and the spectroscopic data: The Sr 4ri states are 
positioned approximately 1.5 eV lower in energy than is 
expected from the experimental spectra as de termined in 
the BIS and XAS spectra (not shown) .122113 The spread 
of the O 2p states is also underestimated by 2 eV. The 
most drastic difference, and one examined many times 
in the literature, occurs in the Ru Ad states, where indi- 
cations of correlations are found. From Fig. [4] it is clear 
that the ti g states at are overemphasized in the LSDA 
calculation. Experimentally the Ru Ad spectrum shows 
only a small weak coherent peak and is 1.5 eV broader 
than the LSDA predicts. As mentioned, the signature of 
strong correlations as observed in the PES data is the 
strong renormalization (or even absence) of the quasi- 
particle peak near e F or large satellite peaks. The large 
coherent t^g peak about 0.5 eV below the Fermi level has 
a substantial incoherent feature near 1.3 eV and is good 
evidence for localized electronic states from strong cor- 
relation effects. We also note that Santi and Jarlborg^ 
suggest that this suppression of the tig states is possibly 
due to small matrix elements for the d — > f and d — > p 
transitions. Using our two "beyond-LSDA" techniques 
to introduce correlation effects, we attempt to address 
whether this incoherent feature, not found in the LSDA 
calculations, is indeed due to strong electron-electron in- 
teractions. 

More recent ultraviolet photoemission spectroscopy 
(UPS) data 5 '- suggests that the Ru stoichiometry plays 
a significant role in determining the spectral weight and 
intensity of the tig peak at the Fermi level. For stoi- 
chiometric SrRuO,3 films, the spectral intensity near the 
Fermi level is in much better agreement with the calcu- 
ated PDOS, while Ru deficient samples have a reduced 
intensity. These facts suggest that previous comparisions 
may have been made with non-stoichiometric Ru sam- 
ples, which may be caused by a high partial pressure of 
oxygen during growth. Although the Ru cation deficiency 
appears to explain the discrepancy with experimentally 
measured spectra, it is difficult to fully remove correla- 
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FIG. 5: (Color online) The total (grey line) and partial spin- 
resolved (shaded) density of states for orthorhombic SrRuC>3 
calculated with U e a = 1 eV (left) and pseudo-SIC (right) are 
shown in each panel, (upper) Sr Ad states, (middle) Ru Ad 
states and (lower) O 2p states. 



tion effects that implicitly result from changes in stoi- 
chiometry, e.g. the d bandwidth can be varied by chang- 
ing the volume of the unit cell via changes in O-Ru-0 
bond angles which occurs upon Ru vacancy formation. 
For example, Siemons and co-workers found the enhance- 
ment of an incoherent peak at approximately 1.5 eV 
below the Fermi level in ruthenium poor samples mea- 
sured with ultraviolet photoemission spectroscopy.^ Fu- 
ture first-principles calculations may be able to identify 
the role of these defects. 

It is clear that the LSDA poorly describes the elec- 
tronic and magnetic structure of bulk SrRuC>3 , and this 
suggests that some underlying physics is missing in the 
local spin density approximation. We next explore two 
extensions of that description to include electron-electron 
correlation effects. 



A. "Beyond LSDA" 

We now explicitly add correlation effects into our elec- 
tronic structure calculations for SrRu03 using the two 
methods outlined above. The pseudo-SIC and LSDA+C7 
methods give very similar results when a U c g — 1.0 eV 
is used. Therefore we first compare the correlated and 
LSDA results, and later point out the small differences 
between results from the two correlated formalisms. 

Orthorhombic Figure [5] shows the densities of states 
calculated with LSDA+t/ with a U e g = 1.0 eV and the 
pseudo-SIC method for Pbnm S1TIUO3. Compared with 
the LSDA, the correlated bands are narrower, with en- 
ergy gaps appearing in both spin channels. The inclusion 
of correlations causes a 70% drop in the total DOS at the 
Fermi level compared with LSDA, with the result that the 
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FIG. 6: (Color online) The orbital bandwidth dependence 
on U e g for orthorhombic SrRuOa. The minority spin states 
are shown by the unshaded symbols and the lines are a 
guide to the eye. Using the pseudo-SIC method, the ma- 
jority (minority) bandwidths are slightly larger than with the 
LSDA+(7 method, but the relative ratios are consistent: Ru 
Ad t 2g = 3.30 (3.50) eV, Ru Ad e g = 4.75 (4.60) eV, and O 
2p = 3.25 (3.50) eV. 



contribution to the Ru 4c? states is almost entirely derived 
from the minority spin channel, and the material is close 
to half-metallicity. This significantly enhances the mag- 
netic properties compared to the LSDA, increasing the 
magnetic moment per formula unit to around 1.0 (J-b, and 
enhancing the exchange splitting of the Ru Ad states at 
T to 0.45 eV («0.30 eV in LSDA). In addition, the peak 
positions of the correlation-included densities of states 
are in better agreement with the experimental spectra 
(Fig. [3|, although the intensity of the O 2p peak at ss 
-7 eV is still too low compared with that of the Ru t^g 
peak. 

The main difference between the electronic structures 
calculated with U c g = 1 eV and the pseudo-SIC meth- 
ods is a larger bandwidth of the occupied orbitals, by 
wl.7 eV, in the pseudo-SIC calculation. This has the 
greatest effect on the Ru Ad bands and the oxygen 2p 
bands near -8.0 eV. As a result, the pseudo-SIC shows 
better agreement with the PES in the bandwidth for the 
O 2p states between -8 and -4 eV, and the ti g states at 
the Fermi level have been more accurately surpressed. 

In order to understand how the hybridization changes 
as electron-electron correlation effects are included, we 
plot in Figure [6] the change in bandwidth for the Ru Ad 
and O 2p states as a function of U c g for orthorhombic 
SrRu03. As the amount of correlation is increased in 
the calculation through the U e g term, the majority Ru 
tig and O 2p bandwidths are strongly reduced, and upon 
narrowing (both by approximately 1.80 eV) half-metallic 
behavior is observed for U c g > 2 eV. On the other hand, 
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FIG. 7: (Color online) LSDA+C/ band structures along T - X 
calculated with VASP for ferromagnetic orthorhombic SrRuOg . 
The horizontal grey line marks the Fermi level. Majority 
(minority) bands are shown as the bold (dashed) lines. The 
highlighted bands indicate the filling of the majority t2 g hole 
pocket, which gives rise to the observed half-metallicity at 
large Hubbard U values. For U c fs = 2 eV, the open circles 
show the results of calculations including spin-orbit coupling. 



only a weak dependence in the orbital bandwidth is ob- 
served for the minority spin states. The valence band- 
width never narrows sufficiently in the bulk material (due 
mostly to the insensitivity of the minority ti g bands to 
correlation effect) to open an insulating gap in both spin 
channels. 

Since the half-metallic ground state that we find 
for [7 ff > 2 eV is not observed experimentally, and 
motivated also by the observation of large magnetic 
anisotropy in Kerr rotation measurements^ on SrRu03 
we have repeated our calculations with spin-orbit cou- 
pling (SOC) effects included. In Figure [7] we plot the 
band structure along r— X in the Brillouin zone as a func- 
tion of increasing U c g. We see that the ti g bands move 
down in energy with increasing f7 e ffj forming a small hole 
pocket which becomes completely filled at U c s = 2 eV 
giving the half-metallic behavior. We note that without 
careful sampling of the Brillouin zone, this hole pocket is 
often missed, and half-metallic behavior can be prema- 
turely predicted. Furthermore, we superimpose the band 
structure calculated with E/ C ff = 2 eV and spin-orbit cou- 
pling in Figure [7| We find here that the degeneracy of the 
tig bands is completely removed, and the highest occu- 
pied majority band is pushed only 0.05 eV higher in en- 
ergy. Furthermore, this splitting decreases at larger U e fi 
values. These results indicate that the inclusion of spin- 
orbit coupling does not have a large effect on the band 
structure. Although spin is strictly not a good quantum 
number, due to quenching of the angular momentum by 
the crystal field, the total angular momentum is well- 
approximated by the spin only component, and therefore 




FIG. 8: (Color online) Calculated spin polarization at the 
Fermi level (Pq F ) for orthorhombic and cubic SrRu03 cacu- 
lated with VASP as a function of U e «. The pseudo-SIC calcu- 
lations give a spin polarization of +2.0% and 8.8% for each 
structure, respectively, (inset) Calculated magnetic moment 
per formula unit for each structure type with LSDA+CA The 
pseudo-SIC calculations yield a magnetic moment of 1.99 and 
1.77/xs/f-u. for each structure, respectively. 

the calculated proximity to half-metallicity in SrRu03 is 
robust to spin-orbit coupling effects. 

Finally we investigate the enhancement of the mag- 
netic properties as C/ e ff is increased. In Figure [8] (inset) 
we show the magnetic moment per formula unit as a func- 
tion of increased correlation for both crystal structures 
of SrRuC>3. For example, the magnetic moment per Ru 
atom is found to be 1.97 [ib with U e g = 1.0 eV and 
1.99 /i b with the pseudo-SIC. At U c g > 2 eV saturation 
of the moment occurs, and the localized ionic moment is 
observed (2.0 hb)- 

Cubic We complete this discussion by describing the 
differences in the bulk cubic electronic structure with cor- 
relation effects added in order to isolate the contribution 
of the octahedral distortions in the orthorhombic struc- 
ture to the bandwidth narrowing. This analysis will pro- 
vide the framework for exploring the Mi-transition in the 
SrRuC>3 thin films. Overall the weight and shape of the 
total density of states with correlations included is con- 
sistent with that calculated in the LSDA, with the excep- 
tion that the large densities of states at the Fermi level 
(majority ti g states) are pushed lower in energy. As was 
observed in the orthorhombic structure, a similar DOS 
is found between both correlation methods, and in gen- 
eral the occupied orbitals with the pseudo-SIC are broad- 
ened by 1.5 eV in energy compared to those calculated 
with the LSDA+[7. As in the orthorhombic structure 
the minority bandwidths in the cubic case are insensi- 
tive to the choice of U e g while the majority O 2p and 
Ru tig bandwidths narrow considerably (1.2 and 1.8 eV, 
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respectively). Furthermore, the observed exchange split- 
tings for the various states are overall larger for the cubic 
structure, and this is consistent with the electronic struc- 
tures calculated with the LSDA. The calculated magnetic 
moment for U c s — 1 eV is 1.64 ps and agrees well with 
that from the pseudo-SIC method (1.77 hb)- Again, for 
values of U c g > 2 eV a half-metallic ground state be- 
comes stable while with the pseudo-SIC a fully metallic 
ground state is always maintained. 

To summarize the bulk SrRu03 results, each of the two 
"beyond LSDA" methods described here improve the de- 
scription of the electronic and magnetic structure. How- 
ever the precise experimental spectra are not fully repro- 
duced although correct peak assignments can be made. 
The addition of a small Hubbard term U e g = 0.6 (1.2) eV 
for the orthorhombic (cubic) structure, or alternatively 
by correcting the SI error in LSDA, also improves the Ru 
tig bandwidths with respect to experiment.^ The total 
width of the O 2p band structure is also increased to ap- 
proximately 7 eV in agreement with the spectral weights. 
We therefore suggest that SrRuC-3 can best be described 
as weakly strongly-correlated. Finally, as stated earlier, 
the intensity at tp has been decreased in comparison to 
LSDA, although it is still larger than experimently ob- 
servered. 



B. Spin Polarization &: Transport Properties 



find that with the planewave code 



-2.95% while 



SrRuOa has been experimentally reported^ to belong 
to the class of negatively spin-polarized materials- char- 
acterized by a greater number of minority spins at the 
Fermi surface which are aligned anti-parallel to the bulk 
magnetization. However, the magnitude of the spin po- 
larization at the Fermi level remains controversial within 
the experimental community, due in part to the differ- 
ent definitions of the spin polarization (resulting from 
the different experimental techniques used to probe this 
quantity), as well as to difficulties in performing the ex- 
periments. Furthermore, the theoretical community has 
also not converged on the magnitude of the spin polar- 
ization, due to the sensitivity of the Ru ti g states near 
tp on the choice of exchange-correlation functional. In 
this section, we perform first-principles transport calcula- 
tions on orthorhombic and cubic SrRu03, and compare 
our results to the available data in the literature. We 
also describe the various definitions of the spin polariza- 
tion commonly used in the literature, and relate them to 
calculated ah initio quantities. 

The spin polarization at the Fermi level Pq f can be 
calculated from the density of states at the Fermi level 
(N tF ) by the following ratio, 







N^ F 



N 



(2) 



Using this definition with the LSDA, the sign of the spin 
polarization for orthorhombic SrRuOs is ambigous: We 



the local orbital code gives a positive spin polarization 
of 2.00%. In constrast, for the cubic structure we find 
a positive spin-polarization (Pq f ) in both first-principles 
calculations, +1.3% (vasp) and +8.8% (siesta). The 
reason for this discrepancy is the sensitivity of the ex- 
change splitting of the Ru Ad bands near the Fermi level 
to the structure. The majority ti g band is positioned 
very close to the band edge, and its precise location is 
sensitive to the finer details of the DFT calculation. As 
a result, the large spread in the calculated spin polar- 
ization is not surprising, and since the magnitudes of the 
spin polarizations are small, changes of a few percent can 
give a change of sign. 

When correlations are introduced, Pg F increases in 
magnitude signifigantly and is negative in all cases. The 
spin polarization as a function of U B s for the orthorhom- 
bic and cubic structures is plotted in Figure [8] The spin- 
polarization calculated for the orthorhombic structure 
with the pseudo-SIC method is -85.7%, compared with 
a value of +2.00% when the SI error is not corrected. 
This value agreesw with that obtained by the LSDA+[/ 
when U e g=l.A eV. We previously found that smaller U c g 
values optimize the agreement between pseudo-SIC and 
LSDA+J7 band structures and magnetic properties, sug- 
gesting that the pseudo-SIC transport results should be 
regarded as providing an upper bound on the magnitudes 
of spin-polarization. For Z7 e ff exceeding a critical value 
of 1.6 eV (3.0 eV) the half-metallic groundstate becomes 
the most stable solution for orthorhombic (cubic) struc- 
ture, as an energy gap opens in the majority ti g band 
and P £F reaches 100%. 

Despite being the most natural definition of spin po- 
larization at the Fermi level, determining Pq f as defined 
in Eq. [2j is a non-trivial experimental process, since the 
spectroscopic measurements required typically have poor 
energy resolution. As knowledge of the degree of spin po- 
larization in a ferromagnet is crucial for its use in spin- 
tronics, several different experimental methods have been 
developed in order to determine this quantity. 

The transport spin polarization can be defined as 



P = 



jt_H 

p + n 



(3) 



where I" is the spin dependent current. However l a 
is not directly observable and must be determined indi- 
rectly. The transport spin polarization now depends on 
the experiment in question, and in particular whether the 
transport is in the ballistic or diffusive regime. In the bal- 
listic limit the current is proportional to N eF vp, while for 
diffusive transport it is proportional to N eF vp (assuming 
both spin species have the same relaxation time), where 
v F are the spin dependent Fermi velocities. 

Therefore the transport spin polarization at the Fermi 
level can be redefined as 



P 



(4) 



10 



Orthorhombic Cubic 



P„ F % (LSDA, pseudo-SIC) 



n = 


+2.00, -85.7 


+8.80, -16.1 


n = 1 


-1.44, -92.9 


-8.99, -50.9 


n = 2 


-15.1, -98.0 


-32.9, -79.5 



TABLE III: Transport spin polarizations calculated with 
SMEAGOL, according to the definition of Eq. fusing both the 
LDSA and pseudo-SIC. Results for both the orthorhombic 
and cubic SrRuOa structures are included. 



where n = 1 fo r ballistic transport or n = 2 for diffusive 
transportPHH. If n — 0, this definition reduces to that of 
the spectroscopic polarization, Pq f . 

An additional definition of polarization is used in 
Meservey-Tedrow style tunneling experiments. Here the 
spin dependent DOS are weighted by their correspond- 
ing tunneling matrix elements. Such an experiment has 
been performed for SRO and report approximately a 
-10% spin polarization.^ Inverse tunnel magnetoresis- 
tance measurements also agree that SRO is negatively 
spin polarizedP This is in agreement with the majority 
of the calculations which find that SRO is a negatively 
spin polarized material at the Fermi surface. 

The point-contact Andreev reflection (PCAR) tech- 
nique, which is based on the process of Andreev 
reflection,^ and developed as an experimental method 
in the work of Soulen et aZP^and Upadhyay et alphas 
been used successfully to determine the magnitude of the 
transport spin polarization, alth ough i t is not sensitive to 
its sign. Experimental results^* 70 * 71 ! using this method 
report values ranging between 51% and 60%. It should 
be noted that in the Andreev experiment the polarization 
is not uniquely defined, in that it must be extracted from 
the data through a fitting procedure and involve terms 
that describe the transmittivity of the interface between 
the ferromagnet and the superconductor. These parame- 
ters are typically difficult to determine precisely and con- 
squently introduce further uncertainty in the experimen- 
tal spin polarization. Furthermore, it is important to 
note that in all PCAR experiments, it is necessary to es- 
tablish whether the transport is in the ballistic, diffusive 
or intermediate regime (non-integer n) which ultimately 
depends on the transmittivity of the interface. The ex- 
perimental results for SRO arc further complicated by the 
fact that the transport in the system has been measured 
in both regimes. 

To allow for a direct comparison with the PCAR exper- 
iments, the transport spin polarization in both the ballis- 
tic and diffusive limit was determined using the ab initio 
electronic transport code SMEAGOLpSlHere we calculated 
the transport at zero bias through both the orthorhom- 
bic and cubic structures and present the results in Ta- 
ble |III| and Figure [9] The shortcomings of the LSDA 
in describing the spin polarization at the Fermi level in 
SrRuOs are again apparent. The highest spin polariza- 
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FIG. 9: (Color online) Spin dependent transport coefficients, 
Nv and Nv 2 , calculated with both the LSDA (unshaded) and 
pseudo-SIC (shaded). 

tion obtained with the LSDA for the orthorhombic struc- 
ture is -15% and it is obtained in the diffusive limit. This 
is notably smaller than the experimental PCAR results 
measuring the same quantity, P^ F . As shown in Fig. [9] on 
changing from Pq f to P| F the polarization increases and 
becomes more negative. Since the group velocity tends to 
zero at the band edge, and is often maximized at the band 
center, higher powers of tl in P^ F suppress the contribu- 
tion of the Ru 4c? states at the band edge while enhancing 
those at the band center. From Figure [9] it is then clear 
that the large negative polarization is a consequence of 
the center of the majority Ru Ad band positioned approx- 
imately 1 eV below the Fermi level, while the minority 
Ru 4ri band center is aligned across the Fermi level. 

Further enhancement is seen by introducing correla- 
tion; for example, by correcting for the SI error, the spin 
polarization increases due to the reduction of the number 
of majority Ru t2 g states at the Fermi level. The corre- 
lated ab initio calculations now give very high spin polar- 
ization, ranging between -85.7% and -98.0% where as the 
highest value achieved experimentally is just 60%. Qual- 
itativity similiar results are found for the cubic struc- 
ture, although the SIC in general has a smaller influence 
on the spin polarization. For example, P{ F goes from 
-8.99% (LSDA) to -50.9% (pseudo-SIC), while P 2 £F goes 
from -32.9% to -79.5%. 

It is also useful to note the strong dependence of spin 
polarization on distance from the Fermi level. In Fig- 
ure [TO] we show for the orthorhombic structure that if 
the Fermi level is moved just 100 meV into the valence 
band, P^ F is decreased to -81.4%, while moving Ep by - 
200 meV decreases it further to -59.8%, bringing it within 
the experimental range of values. In practice this shift in 
the Fermi level can be realized by off-stoichiometric com- 
pounds such as those investigated by Siemons et aZpD 

The discrepancy between the computational and ex- 
perimental results could be then due to a number of fac- 
tors: For example, there are several known limitations 
with PCAR including spin-flip scattering events which 
could drastically reduce the measured value of P tF , as 
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FIG. 10: (Color online) Spin polarization denned according to 
Eq. p] as a function of distance from the Fermi energy (set to 
eV) and calculated with the pseudo-SIC for the orthorhom- 
bic structure. 



well as the possible ambiguous fit of PCAR measure- 
ments to a multiparameter modelP^ We also note that 
spin-orbit coupling, which we did not account for in our 
transport calculations, could reduce the spin polarization 
at the Fermi level. 

Despite these disparities, both ab initio calculations 
and experiment show SrRuC>3 with a high negative spin 
polarization. As expected, LSDA underestimates the 
spin polarization at the Fermi level, whereas the inclu- 
sion of correlation through the correction of the SI error 
with the pseudo-SIC results in much better agreement 
between theory and experiment. 



C. Thin films 

The electronic and magnetic structure of epitaxially 
grown oxide multilayers can be tuned by controlling the 
film thickness. In particular, it has been demonstrated 
that metallic SrRu03 can be transformed into an insu- 
lating state by growing films thinner than five mono- 
layers on SrTiC>3 substrates.^ It was also found that 
the Curie temperature decreases with reduced film thick- 
ness, along with the disappearance of strong ferromag- 
netic order. Photoemission experiments show a shift in 
the spectral weight to the incoherent peak features in 
the spectra, suggesting that these effects are a result of 
changes in electron-electron correlation effects. With our 
first-principles techniques, we systematically investigate 
whether we can reproduce this transition purely from 
structural confinement, or by also including correlation 
effects and/or the octahedral tiltings of the orthorhombic 
structure. For the remainder of this section, we choose 
to include correlation with the LSDA+C/ method, rather 
than the pseudo-SIC method, and note that from the dis- 
cussion so far, both methods reproduce similar electronic 
structures. 

Cubic LSDA slabs To investigate the effects of struc- 
tural confinement on the metal-insulator transition, we 
first performed a series of slab calculations (from 1 to 5 
unit cells thick) on cubic SrRu03 constrained to the cal- 
culated bulk equilibrium SrTi03 lattice parameter. This 
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FIG. 11: (Color online) Magnetic moment dependence on slab 
thickness for cubic SrRuOs. The bulk LSDA magnetic mo- 
ment is shown as the dashed line. The total energy differences 
are calculated for the spin-polarized films with respect to the 
non-magnetic ground state. 



is in part motivated by the fact that good epitaxy is made 
with the substrate surface, and that tilting of the octa- 
hedra may be suppressed. Additionally, it is computa- 
tionally more feasible to systematically investigate these 
smaller supercell slabs. We discuss later the effect of in- 
cluding the octahedra tiltings in the orthorhombic thin 
films; we saw earlier that this structural effect is impor- 
tant in fully describing the subtle details of the electronic 
structure of SRO. In all calculations the slabs were ter- 
minated with a SrO surface, to be consistent with that 
experimentally observed to be the most thermodynami- 
caUy stable.^ 



In Figure 11 we plot the Ru magnetic moment (per 
f.u.) as a function of increasing slab thickness. We find 
that the LSDA films become non-magnetic below a crit- 
ical thickness of only two monolayers; this is lower than 
the experimentally observed loss of the strong ferromag- 
netic order below six monolayers.^ In addition, all of our 
calculations on the cubic films remain metallic down to 
one monolayer. Experimentally the situation is differ- 
ent and insulating behavior is observed in heteroepitax- 
ial thin films at six monolayers on SrTiOs. At one unit 
cell, where all atoms are surface-like, the magnetic mo- 
ment is considerably suppressed from its bulk value, and 
the non-magnetic structure is actually lower in energy. 
An enhancement in the magnetic moment is observed at 
two unit cells in thickness; the moment then decreases 
toward the bulk value as the film thickness grows. For 
films larger than two unit cells, the ferromagnetic ground 
state is always found to be stable. In a mean field the- 
ory approach, the energy difference between the ferro- 
magnetic and the paramagnetic ground states A can be 
expected to be proportional to the Curie temperature T c 
according to ksT c = |A. For the four and five unit cell 
slabs, we find mean field T c 's of 170 and 120 K, respec- 
tively; these values are close to the experimental bulk 
value of 160 K suggesting that even in these thin films 
the strong itinerancy remains. This is consistent with 
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temperatue dependent magnetization data recorded on 
strained and free standing filmiJ^ as well as on ultra- 
thin SrRu0 3 filmsP 1 Additionally, the spin-polarization 
Pq f as a function of slab thickness (not shown) exhibits 
a large negative polarization at two unit cells, while a 
small positive spin-polarization is found with increasing 
thickness (consistent with our bulk spin polarization cal- 
culations). Most importantly, the insulating state is not 
found in any of the cubic slab calculations nor is the non- 
magnetic ground state generally favored (the one unit cell 
case is an exception and is due to competing interactions 
from surface effects). 

To better understand how the magnetism is distributed 
in the slabs, we have also calculated the layer-by-layer lo- 
cal density of states (LDOS). As in the bulk case, on aver- 
age the majority of the spin («65%) is located on the Ru 
atom, with the remaining found on the oxygen network. 
Interestingly, the Ru atoms closest to the surface layers 
experience a suppressed magnetic moment for each slab. 
This is in contrast to most transition metal (non-oxide) 
ferromagnets, where often enhancement occurs due to a 
loss of coordination, weaker interatomic hybridization, 
and enhancement of the orbital angular momentum. In 
this oxide, enhanced covalency at the surface layer may 
be responsible for the reduced magnetism. 

Before adding correlation effects in the cubic slabs, we 
first discuss the changes in the electronic structure due to 
the thin film geometry. The overall shape and weight of 
the density of states for the cubic slab and the bulk cubic 
LSDA calculation are very similar suggesting that con- 
finement effects are minimal. The calculated exchange 
splittings are also similar with the exception that the Ru 
Ad states are split by approximately 0.25 eV. A small 
gap in the minority e g states opens at approximately - 
4.30 eV, and partial occupation of the majority e g states 
occurs; these features are not observed in the bulk cubic 
LSDA calculation. For a free standing, three unit cell 
film we find that the structure has a magnetic moment 
of 1.26 hb and a spin polarization of +13.2% within the 
LSDA, both larger than the bulk cubic values of 1.09 n B 
and +1.3% respectively. The increased positive spin po- 
larization is a result of the band center of the minority Ru 
Ad states shifting to higher energy in the thin films. To 
summarize the results for the cubic SrRuC>3 thin films, 
we do not find a metal-insulator transition as a function 
of film thickness, although we do find a slightly enhanced 
magnetization. We therefore are able to rule out the ef- 
fect of dimensional confinement as the driving force for 
a metal-insulator transition. 

Cubic LSDA+U slabs We now examine the effect of 
adding correlation in the calculations for the cubic thin 
films in order to determine if electron-electron correlation 
in these structures is sufficient to obtain a metal-insulator 
transition. Here we use a U c g — 2 eV, which although 
larger than that we described earlier to more accurately 
reproduce the PES spectra, does allow us to verify that 
in the absence of insulating behavior, the driving force 
for the Mi-transition is not due to intrinsic correlation 



effects. Although the numbers we discuss here are par- 
ticular to a three unit cell thin film we note that the 
general trends are consistent across the series of cubic 
thin films. In contrast to the LSDA calculations, when a 
finite Hubbard U is placed on the Ru Ad states, we find 
that the majority e g states are completely unoccupied, 
and occupation of the majority O 2p states near -2.3 eV 
is enhanced over the minority O 2p states which nearly 
open a gap in the minority spin channael. An enhance- 
ment in the exchange splitting for the Ru d orbitals is also 
observed with U c $ — 2 eV compared to U e s = eV, while 
the valence bandwidth is reduced. For the cubic slab with 
U c ff = 2 eV, the narrowing of the bandwidth nearly sta- 
bilizes a half-metallic ground state, as the majority Ru 
tig bands become completely filled. Despite these small 
changes in the occupation of the Ru Ad levels, we do not 
find an insulating ground state in any of the cubic slabs 
even in the presence of strong correlations (U c s < 6 eV). 
Regarding the magnetic moment in these slabs, we find 

2.0 /is/f.u. for U e s = 2 eV, and a corresponding spin po- 
larization at the Fermi level of -85.9%. These results are 
consistent with the effects of adding correlation in bulk 
cubic SrRuOa, and because an insulating ground state 
is not achieved, we suggest that neither correlations nor 
structural confinement (from our previous discussion) are 
sufficient to induce a metal-insulator transition. 

Orthorhombic LSDA slabs We now address films of 
orthorhombic SrRuC>3 which allow the full octahedral 
distortions found in the bulk experimental structure to 
occur. Earlier we showed that the effect of these dis- 
tortions in the bulk is to reduce the tig valence band- 
width; in this section we examine whether these distor- 
tions with the addition of a confined geometry in a thin 
film form can stabilize the experimentally observed in- 
sulating SrRuC>3 ground state. With the relaxed coor- 
dinates for bulk Pbnm SrRuC>3, we calculate the elec- 
tronic ground state for a three unit cell thick slab sepa- 
rated by 10 A of vacuum on each side and SrO termina- 
tion layers within both the LSDA and LSDA+J7 method 
(U cS = 2 eV)P 

We now discuss the changes in the electronic structure 
of the orthorhombic thin film: In Figure |T2| we show the 
(P)DOS for for the three unit cell slab with and with- 
out correlation. With the LSDA in the thin film system, 
the exchange splittings are similar to the bulk LSDA or- 
thorhombic calculations, and the character around the 
Fermi level remains a mixture of majority and minority 
Ru tig. Energy gaps similar to those found in the bulk 
are observed in other regions of the electronic structure 
with the exception that an additional gap opens in the 
minority ti g states at -2.5 eV, which is not observed in 
the bulk calculation. 

We now compare the magnetic properties of the LSDA 
slab calculation to the bulk Pbnm LSDA calculation. 
With the LSDA method, we find a magnetic moment of 

1.01 fiB per Ru atom and a spin polarization at the Fermi 
level of Pq f — —7.96% (compared to a bulk orthorhom- 
bic structure where a moment of 0.79 ub, and spin- 
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FIG. 12: The total and partial spin-resolved density of states 
for a three unit cell orthorhombic SrRuC>3 slab calculated with 
U e g = eV (shaded) and U c s = 2 eV (unshaded, bold) are 
shown in each panel, (upper) Total (grey) and Sr 4d-states, 
(middle) Ru 4d-states, t2 g and e g , and (lower) O 2p-states. 



polarization of -2.95%). Therefore, we find enhanced 
magnetic properties in the thin film geometry when the 
octahedral tiltings are included. However we still do not 
find an insulating ground state. 

Orthorhombic LSD A +U slabs Finally we incorporate 
correlation into the orthorhombic slab calculations and 
examine the effect on the electronic and magnetic struc- 
ture. We have already demonstrated that a U c h = 2 eV 
is sufficient to establish a half-metallic ground state in 
the bulk orthorhombic structure; therefore, we use this 
limit to establish whether correlation can drive the insu- 
lating ground state. If we do not find a metal-insulator 
transition even at this large Hubbard U value, we can 
be certain that the effect is not due to correlation. In 
general, the shape and weight of the densities of states 
and exchange splittings for the different states remains 
similar to the orthorhombic LSDA slab calculation, how- 
ever unlike the cubic slabs, the valence bandwidth does 
not noticably narrow. 

With the addition of the Hubbard U term in the cal- 
culation, the half-metallic ground state becomes stable, 
with 0.70 eV energy gap opening in the majority spin 
states. This behavior is realized by the minority t2 g 
bands shifting higher in energy while the majority bands 
become completely occupied. The majority e g band is 
also lowered in energy from 1.0 eV in the LSDA slab cal- 
culation (f/ c ff = eV) to 0.50 eV, while the minority spin- 
states move 0.30 eV higher in energy with U c r — 2 eV. 
Similar energy gaps are observed as in the LSDA slab 
calculation, with the caveat that there is no gap in the 
majority O 2p states below the Fermi level; this is due to 



a shift of the O 2p states from the Fermi level to lower 
energy when correlation is added. With a U c g = 2 eV 
we find -100% spin-polarization at the Fermi level and 
a magnetic moment of 2.0 \xb per Ru atom. This effect 
on the magnetism with increased correlation is consistent 
with that found in the bulk calculations. 

In summary, we never find a fully insulating ground 
state in our thin film calculations even in the presence 
of large correlation effects. We have also examined the 
layer-by-layer DOS for each slab (data not shown) and 
have not found an insulating surface layer in any of the 
calculations. However, as a result of the 2D confinement 
in the slabs, we do observe a narrowing of the minor- 
ity t2g bandwidth, and a shift of the Fermi level away 
from the band-center. Furthermore, with correlations 
the Fermi level also is seen to cut across the band-edge. 
These two properties together indicate that SrRuC>3 thin 
films are closer to a metal-insulator instability (with re- 
gards to the bulk), and consequently disorder is more 
likely to induce electron localization and form an insu- 
lating state. From these results we suggest the following 
two possibilities regarding the experimentally observed 
metal-insulator transition: (1) either the transition in 
SrRuC>3 thin films is not an intrinsic property of the 
system, but rather extrinsic and possibly due to surface 
roughness or defects from film deposition combined with 
the band narrowing from confinement and correlation; 
or (2) that SrRuC>3 thin films must be treated with more 
exotic electronic structure methods. 

It is worth mentioning that PES experiments^ of 
SrRuC>3 films grown on SrTiC>3 substrates found a strong 
sensitivity of the t2 g spectral intensity and weight early in 
the deposition process (less than eight monolayers) with 
the intensity of the Ru Ad states becoming strongly en- 
hanced above 15 monolayers. It was found that the film 
growth proceeds with a step terrace mechanism with mi- 
nor atomic diffusion at less than five monolayers and fol- 
lowed by 3D island growth PS The disordered growth pro- 
cess should also reduce the stability of the ferromagnetic 
order, and due to poor percolation pathways, could con- 
tribute to the observed Mi-transition concomitant with 
ferromagnetic ordering at less than five monolayers. The 
disorder at the surface has also recently been compared 
to that at the interface with the substrate (in this case 
SrTiOs) with in situ PES techniques, and it was found 
that the sharp t2 g peak at the Fermi level is greatly sup- 
pressed at the surface, while it persists at the interface.^ 
The decrease in itineracy due to the suppressed DOS at 
the Fermi level was also verified with surface and inter- 
face conductivity experiments. Since our calculations do 
not include any disordered surface configurations or non- 
stoichiometry, future first-principles calculations could 
clarify these competing interactions. We have however 
shown that neither strong correlations nor octahedral 
distortions, nor their combination arc sufficient to repro- 
duce the experimentally observed ultra-thin film metal- 
insulator transition. 
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V. CONCLUSIONS 

Wc have examined the effects of structural distortions 
and correlation effects on the electronic and magnetic 
properties of SrRu03 with first-principles calculations. 
We find that by including weak strong-correlations with 
an effective Hubbard U of 0.6 eV or correction of the 
self-interaction error gives good agreement for bulk or- 
thorhombic SrRuOs with the experimental spectroscopic 
data. The addition of the octahedral distortions leads 
to a narrowing of the majority spin Ru t2 g and O 2p 
states; however the exchange splitting is small with re- 
spect to these bandwidths and consequently a fully in- 
sulating ground state is not obtained. A half-metallic 
ground state was shown to be stable by including moder- 
ate electron-electron correlation effects U > 2 eV, which 
we note has not been observed experimentally. 

The behavior of thin films was also examined in both 
cubic and orthorhombic unsupported films within the 
conventional LSDA approach and with weak correlations 
included. In neither case was the experimentally ob- 
served metal-insulator transition obtained. Since the 
electronic structures of surfaces are very sensitive to 
atomic reconstructions, we suggest that the experimen- 
tally observed metal-insulator transition could be a con- 



sequence of extrinsic defects or an atomically disordered 
surface configuration. 
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